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SUMMARY 


The asymptotic representation is obtained of the atte- 
nuation function in the region of Fresnel's penumbra at ultra- 
short wave diffraction arou.id a planet surrounded by atmosphere. 

The results are presented of numerical calculations for 
an exponential profile of the atmosphere. 

* 

* * 


INTRODUCTION 

The diffraction of electromagnetic waves on an ideally conducting sphere 
situated in a nonuniform medium with dielectric constant profile e(r) = 1 + y 2 /r 2 
was examined in the work [1], 


Analysis of the field in the region of Fresnel penumbra allows us to obtain 
a qualitative representation on the influence of a nonuniform medium on the vari- 
ation of the structure of Fresnel's diffraction field. 

In order to discuss the possibility of applying the described effects in 
case of investigation of planetary atmosphere with emission sources of cosmic 
sonde types' radioeclipse [2], one must be in possession of a solution of an ana- 
logous problem for the profile of dielectric constant, the latter more closely 
corresponding to a certain average authentic planetary atmosphere profile. 

For the latter, the following exponential profile is usually applied: 


( 1 ) 


or 


e(r) = 1 + Aeoe-K'-o) 
n(r) = 1 + A^oe-W'- 0 '; e(r) = ^(r). 


(!’) 


2 


where n(r) is the refraction index. We are interested in a wave band, for 
which the influence of planet's ionosphere may be neglected. 

For the exponential profile e(r) we shall attempt to construct a short- 
wave approximation of the stated problem of diffraction on a sphere in the case 
ka >> 1, where k = 2 tt/X is the wave number, a is the radius of the sphere, and 
to obtain certain numerical estimates. 


1. CONSTRUCTION OF THE SOLUTION FOR THE REGION OF FRESNEL'S PENUMBRA 


In order to construct the solution, let us take the integral expression 
of the radiowave field in a nonuniform atmosphere obtained by Fok [3], and 
applicable for the very general case of arbitrary course of atmosphere's dielec 
trie constant profile from the height e(h) — e(r — a) = e(r) *. 


Then for the Debye potential v of the radial magnetic dioole we shall 
write the expression by means of the contour integral in the form [1] 


v = J-]/ — ‘L— e n^) y v /’(Y, r, r 0 ) f >iv# rfv, 
lb' jtsuiu 


(2) 


where r, r Q are respectively the coordinates of the receiver and of the source; 
■ft* is the angular distance between them, and function F(v, r, r Q ) is expressed 
by the solution of the equation 



(3) 


and satisfies the boundary condition F(v, a, r 0 ) = 0 and the condition of emis- 
sion at infinity. The solution of Eq.(3) may, correspondingly with the idea 
of standard equations method [4] , be represented in the form 


/(v,r) 
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r S(v,r) 
I'i v.r) 




6 '(v.r)- 


«(v,r))[l+0(*a) 

d\/dr | 
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(4) 


where £(v, r ) is the solution of the equation 

kr . i ■— ... 

l K e(p)--(v/p) 2 rfp = v(tgQ-G), cosQsv/g, ( 5 ) 

and the lower limit in the integral krj is determined from the condition 

e(rj ) - (v/kr j ) 2 = 0. 

* It is assumed that in the interval (a, «>) e(r) is not zero and, as r -+• oo 
e(r) 1. 


Applying for Hankel functions H v [ ^ ( r) ] the Debye asymptotic representa- 
tion, and computing the integral in (2) by the stationary phase method, we shall 
obtain for the field an expression valid in the illuminated region and corres- 
ponding to the geometric optics approximation in the solution of the problem. 

For the construction of the solution in the region of Fresnel's penumbra, 
we shall assume for hankel functions the asymptotic expression of the form [5] 

"’ >({)= ~ ■”"*['( f T 1 ''')] (6) 

’ hr, •' 

where Wj (x) is an Airy function; 

y = y(r--« 2 ); 

M= (ka/2)'\ 

The expression for the Debye potential. (2) is then transformed as follows: 


1 i / ~2 ka ( n 2 

IG:iM ’ nr/- 0 siriO l 4 3 

— i y W’} j Vv (0 (»,{v (0 , r , r 0 ]e '^- (t - 
- Y 0 ) { w 2 (t - Y ) - U7, (t - Y ) } dt. 


where 


(p(v, r, r 0 ) = vO + J 4- J — — Jo- ?(v, p) = p 2 e(p) — v* 

Hr, ® fcr, “ 

The complex component _t, over which integration now takes place in (7) , 
is linked with v by the relation [6] 


v = ye(n)/kfl~f — _• Mt., (9) 

Ye («) 

The latter is important when considering the phase funtion <J>(v, r, ro) in 
(8). In the remaining cases one may consider V - ka + Mt. This remark will be 
of use to us in the following. 

For the determination of the form of function G r [v(t), r, r Q ] , entering in- 
to the integral of expression (7) , let us turn at the outset to the physical 
sense of this function. 
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If in (7) we realize the ultimate transition e(p) -*■ 1, we shall obtain the 
expression for the Debye potential corresponding to the case of absence of atmo- 
sphere [7]; postulating G r — 1, i. e. the form of the function is determined 
by the properties of the nonuniform medium surrourding the sphere, namely: it 
describes the refraction attenuation of radiowaves. At the same time, it is dif- 
ficult to determine it from the solution of Eq.(5) in the region of Fresnel 
penumbra. To that effect, we shall make use of the asymptotic solution of the 
diffraction problem on a phase screen (8), valid in the caustics' vicinity. As 
a result, we shall obtain 


C>(v, r, /•„)= -- — }/ 

V«( v - r) q(v, r u ) 


Ji 


(v, r, r 0 ) 


( 10 ) 


where 7? — }'/■* — 2rr 0 cos ft H io 2 , and <{> is differentiated with respect to v(or t) . 
It is not difficult to verify that, as e(p) -*■ 1, G r -*■ 1. 


In the particular case of infinitely remote source (r 0 -*• °°) , we shall ob- 
tain the simpler expressions: 

Gr (v, r, ~ J I V~q(v7r) vV T ( \~r . (11) 

As a slowly varying one, multiplier ✓vG r (v, r, r Q ) may be taken out of the 
integral sign in (7) at t = 0, and we may pass in the expression for the Debye 
potential to attenuation functions introduced by Fok (7], This is fone in detail 
in the work [1]. 


Let us make use of the results of [1] and pass to the particular case of 
interest to us, when one of the correspondents is at infinite range (r 0 -+• °°) . 

Expanding the integrals in (8) in series by powers (v — v^e(a)ka) , we shall 


obtain for 4>(v, r, r Q ) the approximate expressions: 


wnere 


<p(v, r, r 0 ) ez ka r + £r /, 

* r » 

Or — JV«(p)— e(ff)(AVp)*</p 4 - JVc(p) — E(a)(*fl/p) 2 dp; 

hr* hri 

i r ’ hr dp f r * d P W 

^ ^ y^) ( • p*ye(ppV(5)W W* + 1 P : i^j"^)^/T) 2 


( 12 ) 
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It is not difficult to be convinced by means of direct calculation that 
at e(r) — 1 + (y/r)^ we shall obtain the expressions for a r and £ r from [1], 
and for e(r) ■+■ 1, £ r from (13) passes into £ from [8]. 


Passing to the limit for r 0 ■+ 00 , we shall obtain in the illuminated part 
of the Fresnel penumbra the asymptotic expression for the attenuation function 
W [1] in the form 


r 1 sin(Tr 2 — VO ] 


J' 


(14) 


t _ 

where x r == j]Z; r ; n = V Y. 


The expression (14) is fundamental for the discussion of the question of 
the possibility of applying Fresnel difrraction for the study of planetary atmo- 
spheres and conducting numerical estimates. This question is the object of 
discussion in the next section. 


2. INFLUENCE OF PLANETARY ATMOSPHERE ON FRESNEL DIFFRACTION 

As already indicated in [1] and as may be seen from (14), the influence of 
planetary atmosphere is reduced to refraction attenuation of radiowaves, des- 
cribed by the multiplier G r (v, r, r 0 ), and to the variation in the fine struc- 
ture of the Fresnel field of diffraction described by the quantity T r . 

In order to conduct numerical estimates, we shall apply the above obtained 
results to analysis of the diffraction radiowave propagation around planet 
MARS, basing ourselves, first of all, on the fact that an analogous experiment 
already has been realized in the course of the flight of space probe MARINER-4 
[10], Moreover, during propagation in denser atmospheres, such as, for example, 
the Earth's and particularly the Venus' atmospheres, waveguide propagation of 
radiowaves is possible, at which the diffraction phenomenon is getting more com- 
plex. The results of measurements of MARINER-4 have shown that, from the stand- 
point of its electrodynamic properti s, the atmosphere can be characterized by 
the following values of parame ers in (1): Ae 0 = 7.2* 10 — 6 ; 8 — 0.11. The block 
diagram of the experiment and the preliminary computations are described in de- 
tail in the work [2], 

Calculations of the refraction attenuation G r in the atmosphere of Mars 
were conducted by us for several values of Ae 0 = 2An 0 and 8. They are shown 
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in Fig.l by solid lines for certain values of Ae<> — 2An 0 , so as to enable us 
to make comparison with the results of' analogous calculations in the work [2], 
Indicated in Fig.l by dashes and drawn on the basis of radial optics considera- 
tions When computing G r by formula (11), one should note that the differentia- 
tion with respect to V under the integral sign in the expression for <p should 
so conducted that no diverging integral be obtained (to that effect one should 
apply the integration by parts). The calculation of G r was performed by us with 
the aid of a BESM-2M computer. 



Comparison with [2] shows that the 
values of G r obtained in the latter, are 
somewhat greater than our own. This may 
be the consequence of the fact that, gene- 



Fig.l. Refraction attenuation in 
the atmosphere of Mars: 


Fig. 2 


l — fl - 0.05 km- •; i — p - 0,07 K.W- 
3 — fi - 0 , 0 - ' - - - 


1,08 icw-'; 4-0-0,11 *.n 


rally speaking, at propagation of radio- 
waves tangentially to the planet, formulas 
of radial optics are not valid, and this 
is why the calculations made on their basis may lead to somewhat overrated va- 
lues (if one makes in these formulas such approximations that they still allow 
us to conduct calculations at tangential propagation) , which is exactly what we 
observe. This difference may have a significance when interpreting the rv suits 
of measurements of the refraction attenuation for the determination of planetary 
atmosphere's parameters Ae 0 and 0[1O]. 


The results of calculations of G r for Mars' atmosphere parameters, brought 
out in [10] as a function of the dependence of source (or receiver's) raising 
above the planet, h — r — a are compiled in Table 1. 

TABLE 1 

REFRACTION ATTENUATION G r AS A FUNCTION OF HEIGHT h 


‘V *' J J 

—0.810 j 

-0.M5 

— G.vSl 

| 

— !,6tv 

— 3.(£G 

h, KM 

5000 

10 (XHI 

20 OuO 

500.0 

100 1)00 
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Passing to the study of the influence of the atmosphere on the fine struc- 
ture of the diffraction curve, we note that here two effects are possible: on 
account of radiowave refraction there takes place a shadow boundary shift toward 
the side of shadow's region (and alongside with it of the whole Fresnel curve). 

A displacement of oscillations is observed at coincidence with the Fresnel curve 
in case of absence of atmosphere in such a way that their values coincide at the 
boundary of the "geometric" shadow. 

In order to conduct numerical estimates, we shall write 


T r = T-flltt ,-£). (15) 

where T = nS — HM^. \p being the angle of diffraction [9]: 

\p — ^ — a(r) — a(r > ; a * arc cos (a/r) . 

Substituting £ r from (13) into (15), and taking /e(a) - 1 V 2 Ae 0 , we ob- 
tain for the difference £r — C the expression 

Cr - C 2A76u r - Vi W/tfr, (16) 


where ^a r = (a r — a) is the angle of refraction; 


hr 


a r ( r ) l'e(a)ka [ — — — 

* r| P*Ve(p)“ e( fl )(Wp ) 2 


\|), r= 0 — a r (r) — a r (r 0 ), 
or we have for x r the expression 

Tr — t *f- 2 t}.1/5ut — Vi A eo’i 1/vJv. 


(17) 


The position of shadow boundary is determined from the condition T r ** 0, 
\p r being zero at the same time. 


Then, from (17) we obtain the shadow boundary shift in the absence of atmo- 
sphere relative to the position of shadow boundary in the presence of atmosphere: 

t = — 2i vVAcir. (18) 

Formula (18) allows us to determine the angular displacement. The magnitude 
of the latter may be expressed by the distance from boundary position in the ab- 
sence of atmosphere (see Fig. 2: 


1 1 S'. 


( 19 ) 
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Fig. 3. Refraction at 
tangential relative 
to the planet of Mara 
atmosphere radiotrans- 
lucence: 

J — fl- O.OC 2 — p 

- O.OS i — p - 0,0* «M-*; 

4 — [i - 0.07 I — P - 

-0,0 9 * — P - 0,11 km-' 


For the experimental profile of e(r), taking into 
account the parameters characteristics of the atmosphere 
of Mars, the integral, entering into 6a r , may be compu- 
ted approximately. As a result, we obtain 

fior— VsAPoVZ.wjl — ,\ao>2.iap. (20) 

This coincides with [2]. The result of calculation 
by formula (20) is plotted in Fig. 2. 

Besides that, displacement of Fresnel curve's oscil- 
lations are expected at coincidence of the two curves: 
one — in the absence of atmosphere, the other — in the 
presence of it, so that their values coincide at shadow 
bpundary. The magnitude of this displacement -.an be 
determined with the knowledge of "zero" positions 
(points of intersection of oscillations with a straight 
line corresponding to the level |w| — 1 [1]: 


Ti' — tm — 2r\Mba r -{- V?A*onM|>r = t„ — 2i|.WtVi r + '/jAeorpJ/^v. (21) 


Since the quantity i|> r is proportional to the range from shadow boundary, 
we may obtain for the difference V - T nii 6f n ', determining the distance 
between zeros, from Eq.(21): 

fir/ iV„(l + ViApo). (22) 


The relative shift is expressed by the formula 

OF = (Or,/ - 6r„) / 6t m =* VjAeo, (23) 

l.e., the magnitude of the displacement is practically unnotlceable. 


CONCLUSION 

A shortwave approximation for the attenuation function of the field of ra- 
diowaves emitted by a cosmic sonde at: it6 eclipse by the planet is obtained for 
the exponential (or little differing from it) profile of he index of refraction 
of planet's atmosphere, from the point of view of an observer situated on Earth, 
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which allows us t:o estimate the Influence of the planetary atmosphere on the 
Fresnel diffraction of ultrashort radiowaves. 

Estimates were conducted in case of radloecllpses by planet Mars. 

The results of computations show that the described effects, that Is, 
the refraction attenuation, the displacement of shadow's geometric boundary 
and the varlat on of the oscillation period of the Fresnel curve, are small 
as far as the atmosphere of Mars is concerned. Besides, a divergence Is noted 
with the results of refraction attenuation obtained by means of formulas of 
geometric optics. This may have a significance when interpreting the results 
of measurements, the possibility of conducting same is shown in the experiment 
with MARINER-4. 

The author expresses his thanks to N. A. Armand for his constant attention 
in the work and his valuable discussion of the results, to A. M. Ryazantsev and 
V. V. Kirillov for a constructive discussion. 

Manuscript received on 11.7.68 
*** THE END *** 
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